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Abstract 

The product branching ratio, /(b — > Ab) ■ BR(Ab — > AX), where Ab denotes any weakly-decaying 
b-baryon, has been measured using the OPAL detector at LEP. Ab's are selected by the presence of 
energetic A particles in bottom events tagged by the presence of displaced secondary vertices. A fit 
to the momenta of the A particles separates signal from B meson and fragmentation backgrounds. 
The measured product branching ratio is 

/(b -> A b ) • BR(A b -> AX) = (2.67 ± 0.38(stat)t a o f o (sys))%. 

Combined with a previous OPAL measurement, one obtains 

f(b -> A b ) • BR(A b -► AX) = (3.50 ± 0.32(stat) ± 0.35(sys))%. 
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1 Introduction 



In this paper, we present a measurement of the product branching ratio, f(b — > Ab) • BR(Ab — > AX), 
at the Z° resonance. [] In this process a b quark from Z° — > bb decays produces a b-fiavoured baryon 
which decays, directly or indirectly, into a A baryon and other particles. Previous studies of inclusive 
b-baryon decays have emphasized semileptonic decays of the Ab |l|-|3| . 

The result presented here, when combined with the semileptonic branching ratio, fib — > Ab) • 
BR(A b -> A£K), allows one to determine the ratio T(A b -> A£X)/T(A b -> AX) J§. Since T(A b -> 
A£X) depends on V c b and well understood leptonic currents, it may also be possible to extract this 
fundamental weak parameter in a setting with hadronic uncertainties different from those of B meson 
measurements [|],|6|]. 

For this analysis, we select events containing a A particle and a vertex significantly displaced from 
the Z° decay point. This gives a sample enriched in Ab's. Significant backgrounds come from the 
decay of B mesons into A particles, and from b-hadron events where a high-momentum A is produced 
in the primary hadronisation process. These backgrounds are separated from the signal by using a 
simultaneous fit to the momentum and transverse momentum distributions of the A baryon. 

The previously published OPAL measurement of the product branching ratio used a "companion 
baryon technique" to identify jets containing a Ab 0]. That analysis used the momenta of a A and 
an anti-baryon identified in the same hemisphere, whereas this analysis uses the momentum and 
transverse momentum of a A. The two techniques are complementary and have less than 20% of 
events in common. 

After general information about the OPAL detector and Monte Carlo event simulation, we outline 
the event selection which provides an enriched sample of Ab- The backgrounds are then addressed, 
followed by a discussion of signal efficiencies. Systematic errors are discussed in detail. Finally, the 
measured value of the product branching ratio is presented. This is combined with the previous OPAL 
measurement, and is used to update the ratio BR(Ab — > A£X)/BR(Ab — > AX) from Q. 



2 The OPAL Detector and Its Simulation 

The OPAL detector is described in detail elsewhere ]7j . Here, we briefly describe the components which 
are particularly relevant to this analysis. Charged particle tracking is performed by the central tracking 
system which is located in a solenoidal magnetic field of 0.435 T. The central tracking system consists 
of a two-layer silicon micro- vertex detector ||, a high-precision vertex drift chamber, a large- volume 
jet chamber and z-chambers for accurately measuring track coordinates along the beam direction. 

The measurement of specific ionisation in the jet chamber, dE/dx, is used for particle identification. 
Tracks emitted at large angle to the beam direction have up to 159 samplings providing a dE/dx 
resolution of 3.2% §. 

The central detector is surrounded by a lead-glass electromagnetic calorimeter with a wire streamer 
chamber as presampler. The iron magnet yoke is instrumented with layers of streamer tubes which 
serve as a hadron calorimeter and provide information for muon identification. Four layers of planar 
drift chambers surround the hadron calorimeter and serve for tracking muons. 

To obtain momentum distributions of A's from different sources, and to evaluate efficiencies and 
backgrounds, we utilize 6 million Z° — > qq and 3 million Z° — > bb simulated events. The Monte Carlo 
simulation of the OPAL detector is described elsewhere fllQfl. The JETSET 7.4 string fragmentation 



program is used to form hadrons and decay short-lived particles [11,12]. The fragmentation function 



of Peterson et al, is used for heavy flavors, (ep = 0.0038 for b quarks) |13,14[. For this analysis the A 



momentum in the rest frame of B mesons for A's from B meson decay is tuned to match measurements 



by CLEO || 



1 Throughout this paper Ab refers to any weakly-decaying b-baryon. Charge conjugate modes are implied. 
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3 Event Selection 



This study uses a total of 3 554 212 hadronic Z° decays collected by the OPAL detector between 
1991 and 1995. The method for selecting hadronic Z° decays has been described in previous OPAL 
publications [16,17] and has an efficiency of (98.7 ± 0.4)%. 

To select events with a clear two-jet structure, the thrust of the event is required to be at least 



0.8 [18|. Events are also required to be in that region of the detector with good dE/dx and silicon 



micro- vertex coverage by requiring |cos(0y)| < 0.75, where 9t is the polar angle of the thrust axisj^j 
After the angular acceptance cut, thrust cut, and requiring important detector components to be 
operational, N m h = 2 323 302 multihadronic events are retained. 

Reconstructed vertices displaced from the interaction point are used to select Z° — > bb events. 
The primary vertex is determined for each event using the average beam spot position as a constraint 



19, 2CJ. Jets are found using a cone algorithm with a cone having a half-angle of 0.55 radians and a 
minimum jet energy of 5.0 GeV plj| . Both charged tracks and calorimeter clusters not associated with 
a track are used to identify jets. An iterative approach is used when attempting to form a significantly 
displaced vertex, referred to as a secondary vertex, in each jet [22]. Events containing at least one 
reconstructed secondary vertex are retained. The efficiency of tagging a jet associated with a b quark 
is measured to be (21.2 ± 0.9)% for a purity of (95.5 ± 0.5)%. Details of the efficiency calculation are 
described in section [7|. 



4 A Identification 

Events are split into hemispheres using the plane orthogonal to the thrust axis. A particles are identi- 
fied both in the hemisphere with a jet containing a secondary vertex and in the opposite hemisphere 
to increase the sample size. 

The A selection used here is similar to the method described in jl2l| . A's are reconstructed via 
the decay A — > pir. All combinations of well-measured oppositely-charged tracks forming a vertex are 
considered. The higher momentum track is assumed to be the proton. Each track is required to have 
a significant impact parameter in the r-(p plane with respect to the primary vertex to reduce combi- 
natorial backgrounds. The A direction is required to be in the range | cos(#)a| < 0.9. The momentum 
component parallel to the beam line for each track, p z , is re-calculated assuming it originates from 
the reconstructed A decay point 0. A candidates whose invariant mass lies within 8 MeV/c 2 of the 
nominal A mass are accepted if the invariant mass, assuming pion masses for both tracks, does not 
fall within 6 MeV/c of the K° mass. 

The dE/dx identification of the candidate proton requires that the observed energy loss is consis- 
tent with a proton and inconsistent with that of a pion of the same momentum. No dE/dx requirements 
are made on the pion candidate. 

The reconstructed decay point of the A is required to be at a distance greater than 8 cm in the 
r-<p plane from the primary interaction point. There must be no hits in the silicon detector that 
are associated with either of the tracks. The angle 4>(A) between the position vector of the A decay 
vertex and its momentum vector is required to be less than 14 mrad. To reduce A's coming from 
fragmentation, the opening angle of the A direction with the jet axis is required to be less than 0.2 
radians, and the momentum of the A is required to be greater than 5 GeV/c. 

Monte Carlo studies indicate that the fake A's remaining in Z° — > bb events are predominantly 
from real A decays where one decay product of the A is combined with a random track. The fake A 
rate has been studied using side-bands of the A mass distribution and, for the selection criteria of this 
paper, is estimated at 2% in both data and Monte Carlo. The above selection retains 1582 events. 

2 The right-handed OPAL coordinate system is defined such that the origin is at the center of the detector, the z-axis 
follows the electron beam direction and the +y direction points up. The polar angle 6 is defined relative to the +z-axis, 
and the azimuthal angle <f> is defined relative to the +:r-axis. 
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5 Backgrounds 

Besides the Ab — > AX signal, the selected event sample contains the following backgrounds: (1) events 
with a b-hadron and a A produced in the hadronisation process, (2) A's from B meson decay, (3) other 
backgrounds and fake A baryons. These three background sources are discussed below. 

(1) Events where a A baryon arises in hadronisation can be separated from those produced in At, 
decays on a statistical basis. The A baryons from this source generally have low momentum, p, and 
transverse momentum, pt^\ There are two distinct sub-classes to this background: 

• the A is created in associated production with another light baryon within a Z° — > bb event; 

• the A is created in associated production with a primary b-baryon. 

Monte Carlo studies indicate that these two classes of A baryons have similar p and pt distributions, 
which are different from those of the signal. Before the minimum momentum cut, the momentum spec- 
trum for fragmentation A's peaks at 1 GeV, and pt is peaked at 200 MeV. Because these distributions 
have long tails, it is not possible to remove all fragmentation A's with a cut. 



(2) B — > AX decays are a significant background. Although the branching ratio for this process is 
small, the fraction of b quarks which hadronise to B mesons is (89.91^)% in Z° decays Q, so the A 
yield from B mesons is about the same as that from Ab baryons. However, A's from B meson decays 
have a softer pt spectrum than the signal since baryon number conservation requires an additional 
baryon in the decay products of the B meson. 

While not yet observed, it is expected that the Bg meson will produce A's in its decay chain. The 
kinematics are assumed to be similar to B° decays. Excited B mesons are also produced in Z° — > bb 
hadronisation. We assume that these mesons decay either hadronically or electromagnetically to a 
weakly decaying B meson, and that the kinematics are similar to when the B mesons are directly 
produced in the ground state. 



(3) Other backgrounds contribute 3% to the sample. The D + is the only charmed hadron with a 
lifetime long enough to produce a signficant background after requiring a displaced vertex. However, 
since the D + is too light to decay to a A and an anti-baryon, it is only a background when coupled 
with a fragmentation A or a charm baryon to A decay in the opposite hemisphere. The few D + events 
accepted are included in the background class (1) above, since they are kinematically similar. Leading 
A baryons from light quark (u,d,s) decays of the Z° are less than 1% of the sample due to the secondary 
vertex requirement. 



6 Fitting the A Momentum and Transverse Momentum Spectra 

The fraction of A's from Ab — * AX decays is determined by simultaneously fitting the total and trans- 
verse momentum spectra of the A particles. Figure |l] shows the six input Monte Carlo distributions 
used in the fit. They are the p and pt for A's from b-baryons and the two backgrounds: B meson and 
fragmentation. (The uncertainties due to these Monte Carlo distributions are discussed in section ||.) 
The fit returns the fractions of the signal and the two major backgrounds in the data. 

The fit, which takes into account finite statistics, uses a binned maximum likelihood technique 



described in [23]. The 25% correlation between p and pt is not taken into account by this technique. 
The effect of correlations was studied using Monte Carlo samples. This technique underestimates the 
error by 5%, but the central value of the fit fraction is unchanged. We correct the error for this effect. 



3 The transverse momentum of the A is measured with respect to the nearest jet axis. The A is included in the 
calculation of the jet direction. 



5 



Source 


Fit fraction (%) 


Correlations 


A A V 

Ab — > AX 


frag — > AX 


B — > AX 


A A V 

Ab — > Aa 


37.4 ± 5.3 


1.0 


-0.15 


-0.64 


frag — > AX 


37.1 ± 5.1 


-0.15 


1.0 


-0.61 


B -» AX 


25.5 ± 6.6 


-0.64 


-0.61 


1.0 



Table 1: Results of the fit. The fit fractions indicate the fraction of each source in the data sample. 
The correlations between the fit fractions are also shown. 



The fitted distributions are shown in Figure |2| and the fractions are listed in Table [I]. The x 2 of 
the fit is 25 for 27 degrees of freedom. The fraction of signal events is fitted to be (37.4 ±5.3)%, where 
the error is due to finite statistics in data and Monte Carlo. This corresponds to 592 ± 83 Ab — > AX 
signal events. 

The fit method was checked by performing 5000 trial fits on fake data samples which were generated 
by adding the three Monte Carlo sources and allowing the histogram bins to vary according to Poisson 
statistics. The distribution of fit fractions matches the values in the fake data with a standard deviation 
equal to the uncertainty assigned by the fitting routine. 

The stability of the fit result was checked by varying the minimum p and pt allowed in the fit and 
recalculating the product branching ratio for each case. The A momentum was allowed to vary from 
3 GeV/c to 8 GeV/c in 1 GeV/c steps. We also varied the minimum p t cut from 0.0 to 0.75 GeV/c 
in 0.25 GeV/c steps while holding the minimum p cut constant at 5 GeV/c. For all cases the product 
branching ratio was consistent within 10% and there were no deviations that were not compatible with 
statistical fluctuations. 



7 Secondary Vertex Reconstruction Efficiency 

The overall secondary vertex reconstruction efficiency, e b sig , depends on the efficiency of reconstructing a 
secondary vertex in both the unbiased b hemisphere, eb, and the hemisphere containing the decay Ab — > 
AX, ^ax- To measure the efficiency of reconstructing a secondary vertex in an unbiased b-hadron 
hemisphere, we compare the fraction of tagged hemispheres in events with at least one reconstructed 
secondary vertex to the number of events with a reconstructed vertex in both hemispheres. This is 
done for a sample of multihadron events which pass all selection criteria except A identification. 
Solving the following equations for yields the efficiency. 

Rb^h + Rudsctu&sc = flv (1) 

R b el + R udsc e1 dsc = / 2v (2) 
The efficiency of selecting a non-b hemisphere is represented, by the symbol 6 u d SC; 

/i v is the fraction 

of hemispheres with a reconstructed secondary vertex, is the fraction of events in which both 



hemispheres have a secondary vertex, Rf, is r b b/Th a d = 0.2169 |14| and R u dsc = 1 — Rb- This 
measurement yields a value for eb of (21.2 ± 0.9)%, where the error is from finite statistics in the data 
and Monte Carlo, and year-by-year variations in the detector configuration. The effect of tagging 
efficiency correlations between the hemispheres is negligible for this analysis. 

The presence of a high-momentum A in a hemisphere reduces the efficiency for reconstructing a 
secondary vertex. This is true for all sources of A's in the sample. The proton and pion from A's 
in signal events are unlikely to be included in a secondary vertex, which reduces its probability of 
being reconstructed. The A selection enhances the proportion of Ab's which have a shorter lifetime 
than the average b-hadron. Monte Carlo studies indicate that, when the average b-hadron lifetime is 
adjusted to that of the Ab, the shorter lifetime reduces the secondary vertex reconstruction efficiency 
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by a factor of 0.84 ±0.05. Lastly, if the selected high momentum A is from fragmentation, the primary 
b-hadron will have less momentum than usual. A shorter flight distance decreases the probability of 
reconstructing a displaced vertex. 

To calculate the efficiency for reconstructing secondary vertices in Ab — > AX hemispheres, we 
begin by comparing the number of selected A's in the same hemisphere as a reconstructed secondary 
vertex (same-side) to the number of A's in the opposite hemisphere (opposite-side). If the presence of 
high-momentum A's had no effect on the reconstruction of vertices we would expect to find the same 
number of same-side and opposite-side A's in the sample. Instead, we find that the ratio of same-side 
to opposite-side A's is R( data ) = 0.58 ± 0.03 and i?( MC ) = 0.64 ± 0.01. In the Monte Carlo the ratios 
for specific sources of A's are: R ( ^ C) = 0.54 ± 0.03, sfff^ = 0.62 ± 0.03, R^Ln = °- 77 ± °- 03 - 

j{( M C) mus t fog multiplied by a factor of 0.9 to match R( data ) . Assuming that this factor is the same 
for each source, the corrected ratio for At, is R'j^ ata ') = 0.49 ±0.06, where the full size of the correction 
is included in the uncertainty. The efficiency for tagging hemispheres with the decay Ab — > AX is 



4 b ^AX = • < a * a) = (10.4 ±1.3)% (3) 

where the error includes both statistical and systematic uncertainties. 

An event containing a Ab — > AX decay may be tagged by a reconstructed secondary vertex either 
in the hemisphere containing the A or in the opposite hemisphere. The overall efficiency for identifying 
displaced vertices in these events is therefore: 

4 = d b ->AX + eb(l - d b ^Ax) = (29.4 ± 1.5)%. (4) 



8 A Reconstruction Efficiency 

The efficiency of A reconstruction is determined from Monte Carlo for A's satisfying all selection 
criteria. The Monte Carlo simulates well the kinematic properties like A mass resolution, momentum 
distributions, and the dE/dx response of the detector. The overall efficiency for reconstructing A 
particles from Ab decay, e\, is found to be (11.7 ± 0.6)% for a minimum momentum of 5 GeV/c. The 
error comes from Monte Carlo statistics, the 2% fake rate and the tracking resolution. The sensitivity 
to tracking resolution was studied by varying the Monte Carlo resolutions by ±10% which caused the 
A identification efficiency to change by ±0.4%. 



9 Systematic Uncertainty due to Monte Carlo Distributions 

The simultaneous fit to the momentum spectra of the A candidates requires six Monte Carlo input 
distributions: the p and pt of A's from b-baryon, B meson, and fragmentation sources. A systematic 
error is assigned to account for possible mis-modelling of these distributions. Each distribution is 
checked against a data sample, though the comparisons are limited because it is impossible to obtain 
pure, large data samples for the three sources. 

This section describes how the uncertainties on each of the Monte Carlo distributions are deter- 
mined and how they propagate to the fit fraction for the Ab source, f\.. The results are summarised 
in Table §. 



A's from B mesons 

The Monte Carlo momentum spectra of A c 's and A's coming from B mesons are adjusted to match 
CLEO data [15|. In the B meson rest frame, the CLEO data have large errors for A's with momentum 
less than 0.5 GeV/c. A's with low momentum in the rest frame of the B meson are reweighted within 
a range corresponding to these uncertainties, and the fit is repeated for several reweightings in the 
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bources of Systematic hrrors tor jA b 


negative 
errors 


positive 
errors 


p(A) and p t (A) from B — > AX 


-3.5% 


3.5% 


p(AJ ot fragmentation A s 


-1.3% 


0.8% 


Pi (A) of Fragmentation 


-12.5% 


12.7% 


/ A \ J? A A "V 

p(A) from Ab — > AX 


-2.9% 


4.8% 


(A J trom A)-, — y AA 


1 ft ftO/ 


ly.o/c 


Tracking Uncertainty 


-2.7% 


2.7% 


Total 


-21.5% 


24.4% 



Table 2: Systematic errors contributing to the uncertainty in the measurement of the Ab fit fraction, 
Ah" 

selected range, changing the Ab fit fraction by 7%. Values in the center of the range are used in the 
final fit, and an error of 3.5%, half of the variation observed, is assigned for the uncertainty in the B 
meson p and pt distributions. 

A's from Fragmentation 

A powerful technique for isolating A's coming from Ab baryons requires that a lepton with high 
momentum and transverse momentum be identified in the hemisphere with the A 0,0. For these 
studies lepton refers to only electrons and muons. The correlation between lepton charge and baryon 
number of the A is indicative of its origin: combinations with opposite lepton charge and baryon 
number (right-sign) are used to tag Ab — ► A£X events; wrong-sign combinations yield a high purity 
of fragmentation A's [||, which are used as a control sample to compare the p and pt spectra of 
fragmentation A's in the data and Monte Carlo. 

The lepton identification of is used. The A selection is tuned to maximize the number and 
purity of A's from fragmentation in the wrong-sign sample. Differences with respect to the event 
selection of section |] include removing the secondary vertex requirement and requiring the invariant 
mass of the A-lepton pair to be greater than 2 GeV/c 2 to reject A c and B meson decays. With this 
selection, 266 data events are found with a purity of fragmentation A's of 75%. 

Plots a) and c) of Figure || show the comparison of the data and Monte Carlo wrong-sign distri- 
butions. The means of the p and pt distributions for data and Monte Carlo agree well and are listed 
in Table ||. We calculate an uncertainty in the agreement by adding in quadrature the errors on the 
data and Monte Carlo means. An uncertainty of 2.4% is assigned for the momentum distribution and 
4.2% for the transverse momentum distribution. 

To assess a systematic error on fx h due to the Monte Carlo simulation of p or pt, each entry in a 
distribution is multiplied by a factor which raises or lowers the mean of the distribution. For example, 
to see the effect of the 2.4% uncertainty on the p mean we multiply each entry by 0.976, which shifts 
the mean downward. The fit is then repeated with the shifted distribution to observe the change in 
f\ b . The procedure is repeated with a factor of 1.024 to determine the positive error on f\ h . 

The fit fraction f\ h changes by j^J % when the fragmentation p distribution is varied by its un- 
certainty, and by j^il % when the pt distrubution is varied. These changes are assigned as systematic 
errors on f\ h due to uncertainties in the fragmentation A spectra. 

A's from A b 

The technique described in the previous section is also used to compare p and pt distributions for A's 
from Ab — > A£X decays. Right-sign combinations of A's and leptons result in a sample composed mostly 
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Mdata (GeV/c) 


[i MC (GeV/c) 


Sfj, 




WS p(A) 


7.70 ± 0.15 


7.55 ± 0.09 


2.0% 


2.4% 


WS pt{A) 


0.70 ± 0.03 


0.67 ± 0.01 


4.0/0 


4.27o 


RS-WS p(A) 


8.63 ± 0.28 


8.84 ± 0.10 


2.4% 


3.4% 


RS-WS pt(A) 


0.85 ± 0.04 


0.86 ± 0.01 


1.2% 


5.2% 



Table 3: A summary of the A-lepton studies comparing data and Monte Carlo for right-sign minus 
wrong-sign (RS-WS) and wrong-sign (WS) distributions. RS-WS is used to compare A p and pt for 
A's in At, — > A£X decay. WS is used to compare distributions from fragmentation A's. The means 
of the Monte Carlo and data distributions are listed along with their relative difference 5fi and the 
uncertainty on the relative difference which is from the data and Monte Carlo uncertainties added 
in quadrature. 

of signal A's, with the remainder being A's from fragmentation. To isolate the signal shapes, the wrong- 
sign distributions are subtracted from the right-sign. Fragmentation A's and those from light quark 
events populate the right-sign and wrong-sign equally, so the background subtracted distributions 
represent well the momentum distributions of A's from Ab — > A£X decay. 

The event selection used here is the same as for the previous section, except that the minimum 
A-lepton invariant mass cut is set at 2.2 GeV/c 2 , which reduces contributions from B meson and A c 
decays to about one percent of the total. After subtracting wrong-sign from right-sign, the momentum 
distributions in data contain 289 entries. Thep and pt of the subtracted distributions agree well. Values 
are listed in Table || and the distributions are shown in plots b) and d) of Figure |||. 

Shifting the Ab — > AX Monte Carlo distributions by ±3.4% for p and ±5.2% for pt yields systematic 
uncertainties on f\ h of ~t 2 'g % from the p distribution and ligj^ from p t . 

In addition to these tests, we also varied the shape of the pt distribution to ensure that shifting 
the means was a good measure of how variations affect the fit. For several reweightings which skewed 
the distribution to look more like the semileptonic distribution, it was found that the fit result was 
directly correlated to the mean pt . Variations in the mean are a good measure of the sensitivity of the 
fit to variations in the input Monte Carlo distributions. 

10 Other Systematic Uncertainties 

This section discusses systematic uncertainties that have not already been addressed. They are the 
effect of Ab polarisation p5| , Monte Carlo decay model, and tracking resolution. 

Ab polarisation was not simulated in the Monte Carlo. The presence of Ab polarisation shifts 
slightly the momentum distributions of A's coming from Ab's. Possible differences between Monte 
Carlo and data due to polarisation are, therefore, included in the previously assigned uncertainties 
because all Monte Carlo distributions are compared directly to data samples. 

Three additional factors which govern the primary hadronisation process and decay of the Ab and 
its daughter particles could each affect the A momentum distributions. These are the modelling of 
b quark fragmentation and baryon production and decay in the Monte Carlo. Again, since the 
A p and pt for all sources have been compared directly to data samples, the effects are small and 
any contributing uncertainties are included in errors assigned for the modelling of the momentum 
distributions. 

The effect of tracking resolution is also small for this analysis. Varying the Monte Carlo tracking 
resolutions by ±10%, consistent with the known quality of tracking in the simulation, causes the Ab 
fit fraction, f\ h , to vary by ±0.01. This is included in the uncertainty of f\ h in Table ||. 
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11 Consistency Checks of pt for A's from Ab 



Since the overall systematic error is most sensitive to the A pt from Ab decays, we present three addi- 
tional checks on this distribution. While these tests are not as precise as the A-lepton studies described 
in section [j| they are consistent and provide additional evidence that the data is well modelled for 
Ab — > AX decays and that semi-leptonic decays provide a good measure of the Monte Carlo uncertainty. 

(1) One test of the modelling of A p t from Ab — > AX decays is made using the fit directly. The pt 
distribution is shifted until the x 2 °f the fit exceeds the 90% confidence interval. This occurs for shifts 
of +6% and -13%. These limits suggest that the true mean of the pt distribution is somewhere in this 
range. Hence the possible increase in the mean pt for A's from Ab decay cannot be very much greater 
than the assigned uncertainty of 5.2%. 

(2) This test is similar to the A-lepton studies already described in section H but looks for a lepton 
in the opposite hemisphere from the A. This lepton does not bias the inclusive sample of A's, and 
therefore provides a direct check of A's from Ab — > AX decays. When the b quark does not mix in the 
lepton hemisphere, the correlation between charge and baryon number (with the opposite definition 
of right-sign and wrong-sign) still holds for Ab decays. Fragmentation A events still populate the 
right-sign and wrong-sign samples equally. 

Because an invariant mass cut between the A and lepton is no longer meaningful, there are a 
significant number of B meson decays in the right-sign sample. After subtracting the wrong-sign A 
Pt distributions from the right-sign, the Monte Carlo predicts a sample with 65% Ab events, and the 
rest from B mesons. The average pt in data and Monte Carlo agree to within a statistical precision of 
10%. This large uncertainty is due to the many fragmentation and B meson events in both the right 
and wrong-sign samples. While not as powerful as the studies of same-side semileptonic decays, the 
agreement is evidence of good Monte Carlo simulation of non-leptonic Ab decays. 

(3) Finally, the Monte Carlo predicts a Apt mean 13% higher for Ab — * A£K decays than for Ab — > 
AX decays. Three factors contribute to this difference: a 10% correlation between the pt of the 
lepton and A, the A-lepton invariant mass cut, and a smaller multiplicity in semileptonic Ab decays. 
Understanding this difference offers us some confidence that Ab — > AX decays are well modelled. 

The effect of correlation between the lepton and A pt is investigated by varying the minimum 
lepton pt cut and observing the corresponding shift in the A mean pt. For minimum lepton pt cuts 
between 0. and 1.5 GeV/c with 0. as the reference value, a maximum variation of 5% is seen in the 
A mean pt- The Monte Carlo models well the effect in the data. Similarly, the effect of the A-lepton 
invariant mass cut is investigated by varying its value. Once again the Monte Carlo models the effects 
on the data well, with a 7% observed variation in the A mean pt in response to changes in the invariant 
mass cut. Together, the A-lepton correlation and the invariant mass cut account for more than half 
of the 13% difference in the pt means for A's from Ab — > A£X and Ab — > AX decay. 

The remaining difference in the mean pt values can be accounted for by the different multiplicities 
of the two types of decay. Fundamentally, the differences in A pt between semileptonic and hadronic 
decays of the Ab is governed by differences in decay multiplicity. The semileptonic decays are expected 
to have a lower multiplicity, which will cause the mean pt to be higher .0 This effect is investigated by 
examining the track multiplicity in secondary vertices that are in the same hemisphere as the A. 

In the Monte Carlo, the average multiplicity for vertices in hemispheres with a A-lepton pair is 
lower than for hemispheres with just a selected A by 17 ± 2%. In data, this difference is 10 ± 5%. 
The Monte Carlo adequately models the multiplicity difference, providing further evidence that the 
Pt differences between Ab — > AIX. and Ab — > AX events are well modelled 

4 The A momentum in the lab frame is not as significantly affected by decay multiplicity, ft is dominated by boost 
effects. 
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Uncertainty on 
/(b A b ) • BR(A b -> AX) 


N A 

/A b 

Rb 

Nmh 

£A 

e b - 

sig 

BR(A -> pvr) 


1582 ± 40(siai) 
0.374 ± 0.053(stat)1l^ll(sys) 
0.2169 ±0.0012(sys) 
2 323 302 
U.117 ± 0.006(sys) 
0.294 ±0.015(sys) 
0.639 ± 0.005(sys) 


±0.07(stat) 
0.38(siai)l£p(sys) 
±0.01(sys) 

±0.14(sysJ 
±0.14(sya) 
±0.02(sys) 



Table 4: Quantities needed for evaluating the product branching ratio. 



12 Results 

The product branching ratio can be expressed as follows: 



N A f Ah 



f(b - A b ) • BR(A b - AX) = — — (5) 

2R b N mh eAe° ig BR(A -» pvr) 

where JVa is the number of A candidates in the final sample, and f\ h is the fitted fraction of A's from 



A b 's. Rb is the fraction of hadronic Z — > bb decays [14|, N m h is the number of identified multihadronic 



events in the data, e A is the efficiency of reconstructing A — > p7r decays from A b , e b sig is the overall 
secondary vertex reconstruction efficiency for A b — > AX events, and BR(A — > p7r) is the fraction of A's 
decaying to a proton and pion [O. See Table ^ for the values of these quantities. 
Combining these factors, the measured product branching ratio is 

/(b -» A b ) • BR(A b AX) = (2.67 ± 0.38(stat)l[5:6o( s y s )) % - ( 6 ) 

where the statistical error arises from the finite number of data and Monte Carlo events in the fit and 
the systematic error is dominated by the modelling of the A p t spectrum in A b decays. 

To calculate a new OPAL value of the product branching ratio we combine the result presented 
here with an updated value of the previous OPAL measurement using the more recent value of R^ = 

0. 2169 ± 0.0012. This gives, for the previous product branching ratio, (4.00 ± 0A7(stat) ±0.38(sys))% 

01. The event samples for the two measurements of the product branching ratio have less than 20% 
of events in common. Taking into account statistical and systematic correlations, we find 

/(b -» A b ) • BR(A b -» AX) = (3.50 ± 0.32(stat) ± 0.35(sys))%. (7) 

This agrees with and is of higher precision than the value of (2.2^q|)% measured by the DELPHI 
collaboration [^]. 

The value of /(b — > A b ) has been determined to be (lO.llg '\)% [14] using measurements of recon- 
structed A c £ and E£ events. Using this value, we calculate 

BR(A b -» AX) = (35l^)%. (8) 

Finally, we use the data and method of @ to calculate an improved value for the ratio R A £ = BR(A b — > 
A£X)/BR(A b -> AX) = (8.0 ± 1.2 ± 0.9)%. 
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Figure 1: The Monte Carlo p and pt distributions for the signal and the two main background sources 
of b-hadron events containing A's. The six distributions are used as input to the fit. 
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Figure 2: Distributions of A p and pt in data (points with error bars) and Monte Carlo (histogram). 
The three sources: fragmentation, B meson, and Ab are normalized to their fitted fractions. 
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Figure 3: Plots a) and c) represent wrong-sign distributions from A-lepton studies used to compare 
A's from fragmentation in data and Monte Carlo. Plots b) and d) represent the results of the right-sign 
minus wrong-sign subtraction used to compare A's from Ab's. Points with error bars are OPAL data 
and the histogram is Monte Carlo. 
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